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ABSTRACT 


We present a synthesis of fast radio burst (FRB) morphology (the change in flux as a function of 
time and frequency) as detected in the 400-800 MHz octave by the FRB project on the Canadian 
Hydrogen Intensity Mapping Experiment (CHIME/FRB), using events from the first CHIME/FRB 
catalog. The catalog consists of 61 bursts from 18 repeating sources, plus 474 one-off FRBs, detected 
between 2018 July 25 and 2019 July 2. We identify four observed archetypes of burst morphology 
(“simple broadband,” “simple narrowband,” “temporally complex” and “downward drifting”) and 
describe relevant instrumental biases that are essential for interpreting the observed morphologies. 
Using the catalog properties of the FRBs, we confirm that bursts from repeating sources, on average, 
have larger widths and we show, for the first time, that bursts from repeating sources, on average, are 
narrower in bandwidth. This difference could be due to a beaming or propagation effects, or it could 
be intrinsic to the populations. We discuss potential implications of these morphological differences 
for using FRBs as astrophysical tools. 
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1. INTRODUCTION 


Fast radio bursts (FRBs; Lorimer et al. 2007) are j/s-ms dispersed impulses of radio waves that are observed out to 
cosmological distances (see Cordes & Chatterjee 2019; Petroff et al. 2019, for reviews of the phenomenon). Despite a 
growing number of detected bursts, host galaxy associations, and theoretical efforts, the origins of FRBs remain elusive. 
Progenitor models involving cataclysmic events were disfavored by the initial discovery of repeat bursts from an FRB 
(Spitler et al. 2016) and the establishment of repeating sources of FRBs as a class (CHIME/FRB Collaboration et al. 
2019a,b). Many progenitor models include compact objects, with the FRB emission originating inside or outside of 
the magnetosphere, or in shocks in the circumstellar material! (Platts et al. 2019). The recent discovery of a ~MJy 
radio burst from Galactic magnetar SGR 19354-2154 (CHIME/FRB Collaboration et al. 2020a; Bochenek et al. 2020) 
has shown that magnetars are viable candidates as the sources of at least some fraction of FRBs. 

It remains an open question whether all FRBs repeat. So far, it is only possible to decide if any one FRB source 
is repeating by spending considerable telescope time on follow-up observations (Connor & Petroff 2018). Even after 
many hundreds of hours of nondetections, only upper limits on repeat rates can be placed and the matter is further 
complicated by observed orders of magnitude differences in fluence between the initial discovery and repeat bursts of 
some sources (Kumar et al. 2019). Definite proof of the nonrepeating nature of an FRB will likely only be provided 
when an FRB is observed to be associated with a cataclysmic counterpart, such as a supernova or gamma-ray burst. 
By means of repetition statistics one can attempt to answer whether all sources repeat (Caleb et al. 2019; James 
et al. 2020a,b). An alternative means of classifying an FRB as a one-off event or a repeater burst, even if only 
probabilisitically, would therefore be useful. This will aid the selection of the most promising targets for follow-up 
observations and could help uncover the variety of FRBs. The classification of an FRB could, for example, be done 
by way of its spectrotemporal and polarimetric properties or its (position in a) host galaxy. 

Burst morphology (the change in flux as a function of time and frequency) is a powerful proxy for understanding burst 
emission and propagation. FRBs were initially detected in ~ms time resolution, total intensity data, and morphology 
was often discussed in terms of burst structure in band-averaged time series instead of in dynamic spectra of FRBs.? 
Population studies of FRB morphology have become feasible with the advent of large (uniform) samples of FRBs. 
Moreover, with many instruments now capable of saving complex voltage data following the real time detection of an 
FRB, burst morphology can be studied in unprecedented detail for smaller samples of individual bursts as well: at us 
time resolution and with full polarization information. 

Burst structures on timescales of ~ 30 us have been observed for both one-off events and repeater bursts (Farah 
et al. 2018; Michilli et al. 2018). Nimmo et al. (2021) detect burst timescales spanning three orders of magnitude for 
FRB 20180916B, with the briefest bursts lasting only 3-4 ws. This constrains the size of the emission region to be less 
than about 1 km. Comparing Arecibo detections of bursts from FRB 20121102A with the Parkes sample of one-off 
events revealed an apparent difference in bursts widths (Scholz et al. 2017) that was later confirmed by comparing 
repeater bursts and one-off events discovered by the CHIME/FRB project (CHIME/FRB Collaboration et al. 2019b; 
Fonseca et al. 2020). Furthermore, for FRBs 20121102A and 20180916B, burst widths tend to be narrower at higher 
frequencies (Gajjar et al. 2018; Nimmo et al. 2021). At us time resolution, a polychotomy of classes seems to appear, 
between broad singlecomponent one-off events, narrow and multicomponent one-off events, and repeater bursts that 
exhibit downward-drifting chromatic subbursts (Day et al. 2020). 

FRBs have now been detected from 110 MHz to 8 GHz (Gajjar et al. 2018; Pastor-Marazuela et al. 2020; Pleunis 
et al. 2021). Bright apparent nonrepeaters detected by the Australian Square Kilometre Array Pathfinder (ASKAP) 
show knotty, strongly modulated spectra with power concentrated in patches only a few MHz wide (Shannon et al. 
2018). As an ensemble, 23 ASKAP FRBs have an average spectrum that is fitted well by a power-law model with 
spectral index —1.4 or —1.6 reminiscent of average radio pulsar spectra (Macquart et al. 2019). Many bursts from 
repeaters, on the other hand, show Gaussian-like emission bandwidths of 100-200 MHz in the CHIME band (400- 
800 MHz; e.g., CHIME/FRB Collaboration et al. 2019b), and up to a GHz at 4-8 GHz (Gajjar et al. 2018), at varying 
central frequencies. 

Subbursts that drift downward in frequency? were first observed in FRB 20121102A (Gajjar et al. 2018) and studied 
in great detail by Hessels et al. (2019). The phenomenon was also seen in bursts from the second detected repeating 


1 https://frbtheorycat.org/ provides a living review of FRB models. 


? For example, no dynamic spectrum was presented with the first evidence of a multicomponent burst (Champion et al. 2016); but see Ravi 


(2019). 


3 Colloquially known as the “sad trombone" effect. 
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source of FRBs (CHIME/FRB Collaboration et al. 2019a) and now seems ubiquitous among repeaters (CHIME/FRB 
Collaboration et al. 2019b; Fonseca et al. 2020; Day et al. 2020). It has not been observed in all bursts or for all 
repeating sources, but this might be due to the structures being unresolved with limited time resolution or sensitivity 
(Gourdji et al. 2019). Complex structure has also been observed in FRBs that have not yet been observed to repeat 
(see, e.g., Farah et al. 2018; Cho et al. 2020). Wide burst envelopes of up to tens of ms are often comprised of multiple 
spectrotemporal components that can be relatively wide in time in the “wings” of the burst envelope and sharp and 
narrow at the center of emission. The exact appearance of the structures depends on the choice of metric for dispersion 
optimization, but aligning subbursts assuming all frequencies in one component were emitted at the source at the same 
time leads to the least scatter in dispersion measure (DM; Hessels et al. 2019). 

Drift rates, dv/dt (assumed to be linear over the band in which they occur), are measured to be a few to tens of 
MHz ms~! in the CHIME band for all sources where the drift is apparent. FRBs 20121102A and 20180916B so far are 
the only source for which drift has been measured at higher frequencies, and the drift evolves to about GHz ms~! at 
6.5 GHz (Hessels et al. 2019) with an apparent linear trend amidst large scatter (Josephy et al. 2019; Caleb et al. 2020; 
Pastor-Marazuela et al. 2020). Recent wideband and sensitive observations have shown that subsequent subbursts 
spaced < 50 ms can exhibit order-of-magnitude differences in brightness (Chawla et al. 2020; Caleb et al. 2020). 

In summary, FRBs show a variety of morphologies and there are apparent differences between the appearance of one- 
off events and repeater bursts. Any comparisons done so far, however, are between detections with different telescopes, 
with different selection functions, or with a limited sample size. A more systematic study on a large sample of FRBs 
discovered with the same selection function is now possible using the first CHIME/FRB catalog, and is presented 
here. In 82 we summarize the creation of the first CHIME/FRB catalog, and we discuss relevant biases that affect 
observed burst morphologies. In 83 we categorize the catalog into four observed morphological archetypes of FRBs. 
We compare one-off events and repeater bursts statistically in 84. We present subburst separations in 85. We discuss 
our findings in 86 and conclude in 87. 


2. OBSERVATIONS AND ANALYSIS 


The first CHIME/FRB catalog (CHIME/FRB Collaboration 2021, hereafter “Catalog 1") presents all FRBs seen by 
the experiment between 2018 July 25 and 2019 July 2: 61 bursts from 18 repeating sources, plus 474 one-off events.^ 
Catalog 1 contains exposure and sensitivity metrics for all FRB sources during the time frame of the catalog. Events 
are localized based on beam detection S/N ratios, which typically gives a ~15’ precision (for 68% of Catalog 1), but 
for the highest S/N events the localization can be as good as a few arcminutes. In Catalog 1, we analyze the intensity 
data of the beam in which each FRB was detected with the highest S/N. These data have Stokes I total intensity 
at a 0.98304 ms time resolution and with 16,384 frequency channels over 400 MHz. Here, we are not correcting the 
observed distributions of FRBs for the selection function as in Catalog 1. We do exclude three events that were very 
likely detected in far sidelobes from our analysis (see Catalog 1). 

We model all bursts with the least-squares fitting described in Catalog 1, allowing for multiple components. In 
summary, a burst model is described by one DM and one scattering timescale Tscat with a v~+ frequency dependence. 
A burst can consist of multiple components, each with its own intrinsic width and spectrum. The intrinsic widths are 
defined as the standard deviation o of the Gaussian burst that are subsequently broadened by dispersion smearing 
and convolution with an exponential scattering tail (McKinnon 2014). All burst components are modeled at a higher 
resolution than the data (8x in time, 4x in frequency) to model dispersion smearing, before being averaged to the 
instrumental resolution. Based on simulations, we know that we can conservatively measure intrinsic widths and 
scattering timescales at 600 MHz accurately down to 100 s, under the assumption that bursts are comprised of a 
single component that spans the full receiver bandwidth. For high S/N events we can likely get more precise results, as 
long as the models used are an accurate description of the data. Beyond S/N, the precision to which burst parameters 
can be determined depends on the center frequency and bandwidth of the bursts, and the exact dependence will be 
studied in the future. 

CHIME FRBs that are not detected in far sidelobes” have observed spectra that either, i.) extend over the full 
bandwidth or peak at the edge of the band and seem to be best described by a power law, as is common for radio 


4 The repeater bursts in Catalog 1 have all been published before, in CHIME/FRB Collaboration et al. (2019b); Josephy et al. (2019); 
Fonseca et al. (2020). 

5 Detections in far sidelobes show a number of narrow knots of emission (~10 MHz wide each, separated by ~10 to ~250 MHz) due to the 
chromatic beam response in the far sidelobes, see the dynamic spectra of FRBs 20190125B, 20190202B and 20190210D in Catalog 1 for 
examples. In future work we will discuss how the spacing of these narrow features can be used to localize the bursts. 
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Figure 1. Examples of spectral shapes as parametrized by Equation (1) for different combinations of spectral index y and 
running of the spectral index r. Where y is negative and r = 0 the function is equivalent to a power law peaking at lower 
frequencies (green solid curve). Similarly, positive y and r = 0 lead to a power law peaking at higher frequencies. Positive y with 
negative r approximates a Gaussian function (orange dashed curve). With positive y and positive r the result is a steeply rising 
spectrum at the top part of the band (purple dashed curve) and with negative y and positive r the spectrum has asymmetric 
peaks on either end of the band (yellow dash-dotted curve). The vertical axis is scaled linearly. 


pulsars, or ii.) are band-limited and seem to be best described by a top-hat or Gaussian function. The two extremes 
can be relatively well described by one continuous power-law function if we take the common radio spectral index y 
and add an extra “running” term r in the exponent? 


I(v) = Afv v) t el), (1) 


where l(v) is the intensity at spectral frequency v, A the amplitude and vo the pivotal frequency, taken to be 
400.1953125 MHz, the bottom of the CHIME band. 

In the limit where r — 0 this approaches a power-law function. For positive y and negative r this approximates a 
Gaussian function reasonably well. In the other extremes it gives more unrealistic spectra, with a very steep rise on 
either end of the band or with a lack of emission in the middle of the band. Example spectra for possible combinations 
of y and r are depicted in Figure 1. Note that the fitted spectral indices from using equation (1) cannot be directly 
compared to the spectral indices from fits using a power-law function without a running term (where r — 0). The 
bandwidth and peak frequency can be found directly from evaluating the functional form. Alternatively, the peak 
frequency can be found from completing the square of the exponent in the Gaussian form of Equation (1), which 
results in Wertex = vo exp |-y/ (2r)]. 

On top of the burst modeling, we measure burst durations by boxcar convolution with the dedispersed and frequency- 
averaged time series. The measurement of burst duration makes no assumptions about scattering and encapsulates 
multiple subbursts in the case of multicomponent bursts. The width of the boxcar that leads to highest S/N is the 
duration of the burst. The burst bandwidths are derived from the time-averaged burst models, and they are defined 
at the full width at tenth maximum (FWTM) of the peak and capped at the top and bottom of the CHIME band. 


2.1. Biases 


Instrumental and analysis biases, such as the beam response for off-axis detections and limited time resolution which 
leads to unresolved bursts, lead to known biases in the interpretation of FRB morphology. We list various examples 
of bias next. 


2.1.1. Time resolution 


FRBs are us-ms transients with many of them comprised of multiple subbursts with sub-ms widths (Nimmo et al. 
2021). At the 0.98304 ms time resolution of the CHIME/FRB intensity data, the sub-ms structure of bursts will be 


6 This is analogous to the parametrization of scalar fluctuations in the primordial power spectrum used in cosmology (see, e.g., Planck 
Collaboration et al. 2020). There, the power spectrum Ppr (k) = As(k/kg)" 9) with spectral index ns and its first derivative with respect 
to Ink (the “running” of the spectral index) is n(k) = ns — 1 + (1/2)(dns/d1ln k) In(k/ko). Whereas in the cosmological context the spectral 
index and running of the spectral index relate to the physics of the early Universe, the parametrization of FRB spectra here is as-of-yet 
purely empirical, and has no physical underpinning. 
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unresolved and a multicomponent FRB might be interpreted as having just one component. Scatter-broadening can 
also blend multiple intrinsic burst components (e.g., §4.3 in Day et al. 2020). Moreover, any transient search conducted 
at ~1 ms time resolution will will have a lower sensitivity for sub-ms FRBs (Connor 2019). 


2.1.2. Dispersion measure 


DMs measured from data at lower time resolution using incoherent dedispersion (i.e., the data analyzed here) are less 
precise than DMs measured from high-time resolution complex voltage data, where coherent dedispersion is possible. 
The DMs measured at lower time resolution are, however, unbiased (i.e., they are not systematically lower or higher) 
with respect to the DMs measured at higher time resolution; as long as center-of-channel frequencies are used for 
the dispersion correction and frequency channels are sufficiently narrow such that there is very little DM “curvature” 
(i.e., from v-?) across the channel. In the case of narrow frequency channels, the best S/N will occur when all of 
the smeared channels are symmetrically lined up, and the centers of the smeared profiles are defined by the center-of- 
channel frequencies. CHIME/FRB has 16,384 frequency channels over 400 MHz and is in the narrow-channel regime. 
Dedispersion with respect to the bottom or top frequencies of channels will result in a DM bias (of order 107^ pc cm? 
for top-of-the-channel dedispersion of CHIME/FRB intensity data, as determined from simulations). 

What we measure as the DM of FRBs comprised of downward-drifting subbursts depends on the optimization metric 
used. Dedispersed to higher DMs, subbursts can superimpose and lead to higher frequency-averaged S/Ns (Gajjar 
et al. 2018). When optimized for S/N, DMs from repeaters show greater burst-to-burst variation than when they are 
optimized for structure (i.e., when subbursts are aligned; Hessels et al. 2019; CHIME/FRB Collaboration et al. 2019b). 
When FRBs are modeled with multiple Gaussian components, the DMs tend toward values that maximize structure 
(as long as the Gaussians are reasonable approximations of the subbursts). 

When an FRB is comprised of downward-drifting subbursts that remain unresolved, DMs are biased high 
(CHIME/FRB Collaboration et al. 2020b). The exact bias depends on the drift rate, subburst separation and instru- 
mental resolution (Josephy et al. 2019), but we will not attempt to quantify it here. As an example for CHIME/FRB, 
we take two bursts from FRB 20180916B (CHIME/FRB Collaboration et al. 2020b) that show clear downward- 
drifting in high-time resolution CHIME/FRB baseband data and downsample those data to the 0.98304 ms resolution 
of CHIME/FRB intensity data, before optimizing DMs. At high-time resolution these bursts have structure-optimizing 
DMs of 348.78(2) and 348.82(5) pc cm" ?, respectively, but at downsampled time resolution the structure-optiminizing 
DMs are 349.5(3) and 349.3(1) pc cm ?, where the uncertainties are the standard deviation among trials with different 
downsampling offsets, showing that DMs can be overestimated by 0.5-1 pc cm~? in this small sample. We will inves- 
tigate this bias further in future publications of CHIME/FRB baseband data (see Michilli et al. 2021, for an overview 
of the CHIME/FRB baseband analysis pipeline). 


2.1.3. Observed spectra 


'There is a bias due to the chromatic reduction in sensitivity if a burst is detected away from the center of the formed 
beams. Even though the response of all 1,024 beams in CHIME/FRB as a function of frequency and sky position is 
encapsulated in a beam model, the uncertainty on the sky position of each FRB in Catalog 1 is not precise enough 
to meaningfully correct for the beam response. When a source of repeating FRBs is sky-localized to high enough 
precision, as is the case for FRB 20180916B (Marcote et al. 2020), it is known which bursts were detected in and out 
of the beams. Here, as elsewhere, we define the beam edge as the FWHM at 600 MHz of the beams. Comparing the 
bandwidths of FRB 20180916B bursts in Figure 2, it is clear that as a population they tend to have narrow emission 
bandwidths. Moreover, it shows that the narrow bandwidths associated with bursts from this source are real, as they 
are often detected when the beam response as a function of frequency is flat, in the center of the beams. In Figure 2, 
we converted the one-dimensional distribution into probability density functions for visualization, using kernel density 
estimation for data with heteroscedastic errors’ (Delaigle & Meister 2008) as implemented in the CRAN package 
decon? (Wang & Wang 2011). 


3. FAST RADIO BURST ARCHETYPES 


In the observed sample of FRBs in the 400-800 MHz octave at the 0.98304 ms time resolution of the real-time 
CHIME/FRB system, we identify — by eye — broadly four archetypes of FRB morphology: 


T The errors being heteroscedastic, as opposed to homoscedastic, means that not all data points have the same error associated with them. 


8 https://CRAN.R- project.org/package-decon 
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Figure 2. Top: Relative beam sensitivity across the band to bursts from the localized FRB 20180916B at the time of their 
detection, inside (blue) and outside (gray) the FWHM of the formed beams at 600 MHz based on the best available beam model, 
which includes both the primary beam response and the synthesized beams. The rapid variations are due to a “clamping” step 
in the beam forming (Ng et al. 2017). Sensitivity variations due to the polyphase filterbank in the spectrometer are not included 
here and would appear on a much smaller scale. Bottom: Bandwidth FWTM of CHIME/FRB detections of FRB 20180916B, 
separated into bursts that were detected inside (30; blue) and outside (16; gray) the formed beams show that burst emission is 
supressed outside the formed beams. Bandwidths are taken from CHIME/FRB Collaboration et al. (2019b) and CHIME/FRB 
Collaboration et al. (2020b). Both the values (vertical lines) and kernel density estimations are shown (see §2.1.3). 


1. Broadband simple bursts comprised of one peak that can be reasonably well described by one Gaussian function 
in time, convolved with an exponential scattering tail if scattering is not negligible. Their spectra can be well 
described by a power-law function; 


2. Narrowband simple bursts — with spectra that are more like Gaussians; 


3. Complex bursts comprised of multiple peaks with similar frequency extent, with one of the peaks sometimes 
being a much dimmer precursor or postcursor — they can be broadband or narrowband; 


4. Complex bursts comprised of multiple subbursts that drift downward in frequency as time progresses. 


One example FRB of each archetype, as detected by CHIME/FRB, is shown in Figure 3. In Catalog 1, the archetypes 
describe roughly 30%, 60%, 5% and 5% of the bursts, respectively. 

There are multiple ways in which the instrumental response of the telescope can obscure the morphology of a burst as 
compared to how it arrived at the observatory: as described in §2.1.3, beam effects can suppress emission in parts of the 
band, making a burst appear narrower in frequency (type 1 — 2) and potentially concealing drifting subbursts (type 4 
— 3). Bursts with structure below the time resolution or with low S/N will be unresolved (type 4 — 2). Furthermore, 
in the context of this paper, “broadband” is used to describe emission that extends roughly from 800 MHz down 
to 400MHz; however, it does not necessarily imply that the burst emission extends over many GHz of bandwidth. 
Moreover, drift rates => 400 MHz ms~! (receiver bandwidth over time resolution) will not be identified as drifting, 
appearing as type 3 instead of 4 (see also the discussion on drift rate measurements in CHIME/FRB Collaboration 
et al. 2019b). 

A broadband and singlepeaked morphology has so far never been clearly observed in a repeater burst detected by 
CHIME/FRB (see dynamic spectra in CHIME/FRB Collaboration et al. 2019a,b; Josephy et al. 2019; Fonseca et al. 
2020; CHIME/FRB Collaboration et al. 2020b) and downward-drifting subbursts seem ubiquitous among repeater 
bursts (detected in at least 19 bursts from 13 different sources among the 58 bursts from 19 sources in CHIME/FRB 
Collaboration et al. 2019a,b; Josephy et al. 2019; Fonseca et al. 2020). 
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Figure 3. Dedispersed dynamic spectra of the four FRB morphology archetypes described in the text, as detected by 
CHIME/FRB. These are FRBs 20190527C, 20190515D, 20181117B and the 2019 August 10 burst of repeating FRB 20190117A. 
Knots of intensity in the burst spectra are instrumental in origin. 


4. REPEATERS VERSUS NONREPEATERS 


In Figure 4, we compare the morphology of FRBs that have not yet repeated with those of FRBs from repeater 
sources using the fitted intrinsic widths and spectral runnings as a proxy for bandwidth. We consider the sample of all 
repeater bursts as well as the sample with only the first repeater detections, as the former is biased by the detection 
threshold for subsequent bursts from the same source being lower and by containing a disproportionate number of 
bursts from more prolific repeater sources. In Figure 5, we compare the durations and bandwidths of FRBs. This is a 
slightly more intuitive comparison as the two values can be more easily estimated by eye from a dynamic spectrum. 

By eye the parameters associated with one-off events and repeater bursts already appear to be distributed differently. 
Even though burst spectra are described by two parameters (the spectral index y and running r, or bandwidth and peak 
frequency), we compare the parameter distributions separately, as nonparametric statistical tests that ask whether two 
samples are drawn from the same probability distribution are ill-defined in more than one dimension (e.g., Feigelson 
& Babu 2012). 

We take the intrinsic widths, burst durations, spectral runnings and burst bandwidths and compare all one-off FRBs 
with the first detected bursts from repeaters, using the k-sample Anderson-Darling (AD) test (Scholz & Stephens 
1987), as implemented in the CRAN package kSamples.? We do not include a comparison of the spectral indices as 
this parameter is correlated with the spectral running and the running more directly distinguishes broadband (r ~ 0) 
and narrowband (r « 0) emission. To account for the uncertainties on the measurements of the intrinsic width 
and spectral running we Monte Carlo (MC) resample 10,000 times, each iteration drawing a measurement randomly 
from its Gaussian uncertainty region before recalculating the p-value from the AD test comparison. Across all MC 
iterations we calculate a 95% confidence interval upper limit on the p-values. We repeat the analysis after imposing 
a detection S/N > 12 threshold for inclusion (below this threshold real events that were detected in the search might 
have been misclassified as noise or radio frequency interference by humans verifying candidate FRBs). The results of 
the analyses, compiled in Table 1, show that among all four variables there exists a statistically significant observed 
difference between the two samples. It is thus clear that one-off events and repeater bursts differ strongly in average 
width and bandwidth. 


9 https://CRAN.R-project.org/package=kSamples 
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Figure 4. Spectral runnings and intrinsic widths from model fits to the FRBs in Catalog 1, separated in one-off events (green 
diamonds), repeater bursts (orange open circles) and first repeater detections (purple open circles). Error bars represent lo 
uncertainties on the least squares model fits. On the top and on the right the 1D distributions (vertical/horizontal lines) and 
associated kernel density estimates (see §2.1.3) of the respective parameters. 


5. SUBBURST SEPARATIONS 


For FRBs that have multiple components, the difficulty in asserting a single event (under a single “envelope” ), as 
opposed to multiple events close in time, grows larger as the subburst separation increases. Using the model fits from 
Catalog 1, we can test whether there is a bimodality in subburst separations (i.e., if the data support a threshold 
separation below which two subbursts are part of the same FRB envelope and above which two subbursts are instead 
two separate FRBs). At the same time, we can also investigate the separation between the frequency center of 
subsequent components. 

We calculate the time separation At = TOA;,; — TOA; and the frequency separation Av, = vci41 — vc; where 
TOA, is the time-of-arrival and ve, the frequency center of component i. A negative Av, implies a downward drift. 
Note that there is an unquantified uncertainty in the peak frequency, especially for those subbursts that extend beyond 
the edges of the receiver bandwidth. As a separation is only defined for FRBs with two or more components, we only 
include those FRBs in this analysis. An FRB comprised of n subbursts results in n — 1 data points in this analysis. 
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Figure 5. Bandwidths and durations of the FRBs in Catalog 1 with normalized histograms on the sides. FRBs are separated 
into one-off events (green diamonds) and repeater bursts (orange open circles). The panels on the right show smoothed and 
normalized distributions of all one-off events (green) and all repeater bursts (orange), respectively. Only bursts with detection 
S/N > 12 are included. Note that the one-off events may contain as-of-yet undiscovered repeaters sources (e.g., due to limited 


exposure or source activity). Note also that a subset of bursts have underestimated bandwidths due to detections away from 
the center of a formed beam (cf. 82.1.3). 


Table 1. One-dimensional statistical comparisons of one-off FRBs and repeater bursts, including all FRB subbursts (*individ- 
ual") or per-FRB averages of subburst parameters (“averaged”). Comparisons are done without and with a cutoff on S/N = 
12, with the results for the samples with a S/N cutoff added in parentheses. 


One-off events vs. Subbursts Sample size p-value 


Width? Duration Spectral running? ^ Bandwidth 

(6) 10? (1073) 1077 (1078) 107? (1072) 1079 (1075) 

(0) 1074 (107°) 1074 (1074) 107° (1077) 107? (107°) 

vs. 93 (57) 10-P (107'') 107?! (107?) 10786 (107°) 107” (1072) 
(33) 107+? (107°) 10-!? (10779) 10-74 (10716) 1071? (1075) 


Repeater firsts Individual 503 ( 
Repeater firsts Averaged 471 ( 
Repeater bursts Individual 503 (342 
Repeater bursts Averaged 471 ( 


495% confidence upper limit from 10,000 Monte Carlo trials. 


In Figure 6, we compare all subburst separations. From the distribution of subburst separations in time there is 
no clear bimodality between short (here ~ms) and long (here ~tens of ms) separations in this sample, with only two 
bursts with separations > 30ms. The data seem consistent with being drawn from an exponential-like distribution 
with a long tail. This holds true for both one-off events and repeater bursts, with repeater bursts on average seeming 
to have larger subburst separations. However, note that the ~1 ms time resolution is left-censoring the distribution. 
Considering the frequency separations, all repeater bursts have Av, < 0, which means that they do not drift in 
frequency or that they drift downward in frequency.!? 

Some of the one-off events might turn out to be bursts from as-of-yet undiscovered repeater sources. The overlap 
of one-off events with the sample of repeater bursts may be due to undiscovered repeaters.!! Three one-off events 
have apparent upward drifts between some of their subbursts, but all three are also broadband bursts, where the 
determination of the center frequency within the band is ambiguous and the uncertainty on the center is difficult to 
quantify. The frequency separations for those bursts are thus likely consistent with being zero. The dynamic spectra 
of these FRBs, 20181226B, 20190122C and 20190131D, are shown in Figure 7 for reference. 

For comparison, we have also added the CHIME/FRB-detected burst from SGR 19354-2154 to Figure 6. The two 
observed components for this burst clearly have different spectra, with the first peaking at the bottom of the CHIME 
band and the second at the top. This results in an upward drift of the frequency centers. The two components are 
separated by ~29 ms and could be independent bursts rather than subbursts of the same envelope (CHIME/FRB 


10 Note though that bursts 24 and 25 from FRB 20180916B (not included in Catalog 1) are separated in time by only about 60 ms and appear 


to have an upward drift of about 100 MHz (CHIME/FRB Collaboration et al. 2020b). 
11 We select the best candidates based on burst morphology and exposure/burst rate in Appendix A. 
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Figure 6. Subburst separations in time and frequency in Catalog 1 with bursts from apparent nonrepeaters (green diamonds) 
and repeater bursts (orange and yellow open circles). Note that only FRBs that have two or more subbursts are included. Of 
the 18 repeater bursts, 8 are from FRB 20180916B. Those 8 have been colored yellow in the central panel to avoid biasing the 
interpretation. The CHIME/FRB detection of a burst from SGR 1935+2154 (purple star) is shown for reference. On the top 
and on the right the 1D histograms of the respective parameters, including all repeater burst pairs in orange. 
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Figure 7. Dynamic spectra of FRBs 20181226B, 20190122C and 20190131D with three, six and two subbursts, respectively. 
These three FRBs show apparent upward drifts between subbursts in Figure 6 (see 85). The band-averaged time series and 
time-averaged spectra are shown in the panels on the top and on the right, respectively, with the fitted models overplotted in 
a blue solid line. The burst durations and frequency extent are indicated by a blue shaded region. The burst names and DMs 
(in pc cm? units) are displayed in the panels with the time series. Data are from Catalog 1. 

Collaboration et al. 2020a). This burst thus seems to be an outlier among the FRBs in Catalog 1 and its morphology 
is not standard for repeating FRBs. It should be noted, however, that this magnetar burst had exceptionally large 
observed flux density due to its proximity and that for a lot of more distant FRBs a second (or n") component might 


be buried in the noise. 


6. DISCUSSION 


We find an observed difference in the burst morphologies of one-off FRBs and repeater bursts. It is unlikely that 
these differences are caused by instrumental effects, such as beam response towards the sky location of an FRB at 
the time of its detection or the time resolution of the data, because both populations are subject to nearly identical 
biases. Similarly, the CHIME/FRB selection function that Catalog 1 presents, based on injections of simulated events 
in the real-time system, affects both one-off event and repeater burst detections. The selection function shows that 
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CHIME/FRB is only recovering a fraction of wide and highly scattered events, which implies that the differences in 
width/duration between the two populations might even be understated in the current sample. 

Moreover, it is clear that many FRBs are not simply broadband, narrow impulses that are potentially scatter- 
broadened but that FRBs are often comprised of multiple subbursts. This becomes especially apparent at sub-ms time 
resolution (e.g., Farah et al. 2018; Day et al. 2020). It can be difficult to disentangle intrinsic burst structure and 
propagation effects, especially for FRBs detected with limited time resolution or S/N (Gourdji et al. 2019). This in 
turn, hinders the precise determination of DMs, scattering timescales and subburst separations. 

Next we discuss the potential interpretations and implications of burst morphology for FRB models and applications, 
respectively. 


6.1. Potential interpretations 


It is possible that a selection effect due to beamed emission causes the observed difference in bursts widths, assuming 
there exists a positive correlation between burst width and beaming angle (Connor et al. 2020). A simulation by Connor 
et al. (2020) of this putative effect reproduces the difference in burst widths between one-off events and repeater bursts 
qualitatively. The model predicts a correlation between repeat rate and burst width that cannot be significantly 
detected with the currently published set of repeater sources. An explanation for the narrow bandwidths of repeater 
bursts in this framework is currently lacking, and not accounted for in any model to our knowledge. More generally 
(i.e., without necessarily invoking beaming), the two populations could be manifestations of two different regimes out 
of a continuum of source activity that is correlated with the burst morphology. 

The differences might also be due to propagation effects. Cordes et al. (2017) proposed plasma lensing as a way to 
explain the complex time/frequency patterns, which we now recognize as typical of the repeating FRB population, in 
the first discovered repeater, FRB 20121102A (e.g., Hessels et al. 2019). Cordes et al. (2017) showed that impulsive 
radio waves passing through discrete AU-scale structures can experience strong frequency-dependent focusing that 
results in high-gain caustics that alter light curves and spectra. They suggested such structures could arise from 
one of three sources: a source-induced bow shock in a dense medium, filaments in a surrounding supernova remnant 
(SNR), or filaments in a surrounding H II region. All these options tend to implicate environments like those in the 
disks of spiral galaxies or star-forming regions. However, if this were the correct interpretation for the time/frequency 
structure of repeaters, then it suggests apparent nonrepeaters must not be as commonly subject to plasma lensing, 
which argues for them having different environments. Yet, as shown in Catalog 1, repeaters and apparent nonrepeaters 
have DM distributions that are consistent with arising from the same underlying distribution. The same is true of their 
scattering time distributions. This suggests similar environments, at odds with the plasma lensing implications, unless, 
for example, the two classes represent only an age disparity, with repeaters somewhat younger (hence preferentially 
inside SNRs), but having largely similar galactic environments. Though there are small disparities between the host 
galaxies of the two classes, robust conclusions regarding their differences cannot yet be drawn due to the small number 
of repeating sources with host galaxy associations (Bhandari et al. 2020; Heintz et al. 2020). We note that the strong 
preference for downward-drifting subbursts in repeater bursts is problematic in the plasma lensing model as already 
noted by Hessels et al. (2019); our results only exacerbate the problem. 

A different possibility is that if repeater sources are preferentially located in binary systems (cf. the periodic activity 
of FRBs 20180916B and 20121102A; CHIME/FRB Collaboration et al. 2020b; Rajwade et al. 2020; Cruces et al. 2021) 
the intrabinary or circumbinary material might cause a propagation effect that leads to the wider bursts and narrower 
bandwidths. If so, then periodicities in repeat bursts should eventually be shown to be common. Plasma lensing has 
been observed in Galactic binary pulsars (Main et al. 2018; Bilous et al. 2019). 

Finally, it is possible that the observed differences are intrinsic, with one-off events and repeaters forming two 
distinct populations of FRBs produced by different progenitors and/or emission mechanisms. However, again, it 
would be puzzling that no differences have so far been seen between the observed DM and scattering distributions 
of the two populations if they have different progenitors (Fonseca et al. 2020; Catalog 1). Perhaps one-off events 
and repeaters are produced through the same astrophysical formation channels, with the differences in burst rate and 
complexity of the FRB emission being due to a separation in age, magnetic field (strength and/or orientation) or some 
other parameter intrinsic to the source. 

Some proposed models for FRBs predict narrowband emission: decelerating blast waves result in Av/v ~ 1 emission 
(Metzger et al. 2019), which is narrowband, but still necessitates some external effect to produce the Av/v ~ 0.1 that 
is observed for repeater bursts. Superradiance in a molecular cloud triggered by a coherent pulse from, e.g., a pulsar, 
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appears to be able to explain repeater burst bandwidths and spectral structure (Houde et al. 2019). These models, 
however, do not predict a difference between one-off events and repeater bursts. The interaction of hypothetical axion 
stars with neutron stars and the accretion disks of black holes has been proposed to lead to one-off events and repeating 
FRBs, respectively, with different center frequencies and bandwidths of emission between the two populations (Iwazaki 
2020). 

Repeater bursts might be made up of multiple subbursts analogous to the main pulse of the Crab pulsar (see, e.g., 
Cordes & Wasserman 2016), with many of the subburst separations below the time resolution of the CHIME/FRB 
~1 ms data and most other detection systems. If some physical process causes a distribution of subburst separations 
that peaks below the time resolution of an instrument, it is most likely to detect burst separations close to the time 
resolution of that instrument, as can be observed in Figure 6. Mapping out the precise wait-time distribution will help 
differentiate between FRB emission models and, e.g., the CHIME/FRB baseband data will help address this. 


6.2. FRB classification 


Our results show that FRBs can be probabilistically classified as either a one-off event or a repeater burst, based 
on only their burst morphologies. In Appendix A we present one-off events from Catalog 1 that are likely to yield 
repeat bursts, based on their morphology. As a next step, a predictive model can be built using a machine learning 
algorithm, to investigate how accurately this classification can be performed with the current data set. These kind 
of classifications can be made more precise in the future by not only considering the burst morphologies, but also, 
e.g., the brightnesses, distances, polarization information and host galaxy associations of the FRBs. A multiparameter 
comparison will also likely reveal if there exist more than two observed classes of FRBs. Concurrently, it might be 
possible to train a classifier on the FRB dynamic spectra directly, as has been successfully done to unambiguously 
separate short and long GRBs (Jespersen et al. 2020). From expanding the CHIME/FRB sample of FRBs and a 
more quantitative understanding of the experiment’s selection effects we will soon be able to determine if the observed 
classes of FRBs are distinct or instead various realizations of one continuum of FRB properties. 


6.3. Implications for FRBs as astrophysical tools 


FRBs have been proposed as probes for solving long-standing open questions in astrophysics, such as weighing the 
“missing” baryons in the intergalactic medium (e.g., McQuinn 2014; Macquart et al. 2020), measuring the reionization 
history of the Universe (e.g., Inoue 2004) and the optical depth of the cosmic microwave background (e.g., Fialkov 
& Loeb 2016). These applications depend on measuring FRB DMs and redshifts. Even though we have shown that 
burst morphology limits the determination of DMs to ~0.1 pc cm"? precision in many cases"?, this will not affect the 
use of FRBs for these applications, as the error budgets for these analyses are dominated by the uncertainty on host 
galaxy contributions to DMrng that are of order 100 pc cm-?. 

'The progenitors and local environments of repeating FRBs can be probed by DM evolution of the bursts. FRBs are 
often theorized to live in dense environments, like pulsar wind nebulae or supernova remnants, especially in models 
in which the FRB source is a young neutron star. In those cases, the rapidly changing environment can introduce 
significant DM changes up to a few tens of pc cm ^? over the first few years after the supernova (Piro & Gaensler 
2018), that can be probed with ~0.1 pc cm"? precision. It is also possible to probe the kind of orbital DM variations 
due to changes in the line-of-sight electron density observed for binary Galactic pulsars with B (up to 4 pc cm?) or 
Be (up to 11 pc cm?) type companions (Madsen et al. 2012; Wang et al. 2004) but not the 107? pc cm? variations 
for pulsars with brown dwarf companions (Main et al. 2018). It has furthermore been calculated that the propagation 
of luminous FRBs (isotropic energy of 10*° erg s71) can accelerate electrons in the local environment (< 1 pc) of an 
FRB source, which can lead to burst-to-burst DM variations up to 107? pc cm~? depending on the wave strength of 
each burst (Lu & Phinney 2020). These kind of measurements can only be made with better DM precision at ~ ps 
time resolution (Cho et al. 2020). 

FRBs can potentially be observably gravitationally lensed and probe dark matter distributions in the Universe (e.g., 
Munoz et al. 2016; Sammons et al. 2020). Whether wave effects are directly visible in the burst morphology of an 
FRB depends on the Fresnel scale of the lensing scenario, which in turn depends on the physical parameters of the 
lens and the observing frequency. In the infinite frequency (“geometric optics") limit (or if FRB images decohere due 
to the presence of a scattering screen in the FRB’s line-of-sight), coherent effects due to the wave nature of light are 


1? Though note that the achievable precision depends on the observing frequency and the resolution and sensitivity of the instrument used. 
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unimportant. In this case two lensed copies will appear as two bursts with similar morphologies and polarimetric 
properties, but with different magnifications. In the finite high frequency (“eikonal”) limit, interference patterns 
typical of a double-slit experiment may be seen. In the low frequency (“wave optics”) limit, effects of diffraction may 
be observable as well (Jow et al. 2020; Feldbrugge & Turok 2020). 

A first hint that two FRBs are lensed copies of each other could thus be similar burst morphology and sky location. 
An unambiguous detection of lensing can be made when the electric fields of two FRBs correlate. In the geometric 
optics limit, the detection of a lensed event could be claimed based on solely the similar burst morphology and 
polarization (on top of similar sky location), in the absence of strongly chromatic effects in the vicinity of the lens. 
The morphological diversity among especially repeater bursts could make the identification of potentially lensed events 
easier (i.e., it would be more difficult to assert two events are the same if all events were identical). 


7. CONCLUSIONS 


We have considered the burst morphologies of all FRBs in the first CHIME/FRB catalog using simple model fitting 
to their dynamic spectra. We have compared the widths and bandwidths of one-off FRBs with those of repeater bursts. 
There is a statistically significant observed difference between the two, with repeater bursts on average having a longer 
duration and being narrower in bandwidth. 

This difference could be due to a beaming or propagation effect, or it can be intrinsic to the populations. We have 
also considered subburst separations in the first CHIME/FRB catalog. We find no clear bimodality in separation 
between subbursts that clearly fall under one envelope and those without a “bridge” in emission between subbursts in 
this sample. 

Looking forward to a second CHIME/FRB catalog, better burst-fitting techniques will lead to more accurate model 
fits and thus a better morphological population comparison. For example, the beam model and sky position uncertainty 
regions could be taken into account to get a better handle on the bandpass response and hence the burst bandwidths. 
FRB subburst separations can be investigated further by using CHIME/FRB baseband data that are saved to disk 
for about a hundred FRBs that are part of the first CHIME/FRB catalog. These data, with a time resolution of 
2.56 us, can resolve a lot of the structure in FRBs that remained unresolved in this analysis. It will potentially reveal 
a plethora of burst structures below a ms. Additionally, the baseband data can be phased up towards the best known 
sky position of an FRB, which results in data that are unaffected by the synthesized beam (though not the telescope 
primary beam) that causes much of the instrumental response that affects the intensity data (Michilli et al. 2021). 

Many one-off FRBs span the full CHIME band, but it is not clear how broadband their spectra really are. Future 
very-wide-band detections of one-off FRBs or detections with multiple facilities at the same time will reveal how far 
the bursts’ instantaneous spectra reach. 

Finally, we see two possible avenues for solving the conundrum of the significant difference in repeater versus appar- 
ent nonrepeater burst properties, in spite of the two classes’ similarity in DM and scattering time distributions. One 
may be via comparative studies of low-DM sources, which, due to their proximity, could eventually reveal differences 
in multi-wavelength properties which could solidify distinct natures. Alternatively, a large number of sub-arcsecond 
FRB localizations could eventually help pinpoint locations within host galaxies and show differences between repeat- 
ing sources and one-off events that will support different source classes. This should be enabled by the upcoming 
construction of CHIME/FRB Outrigger telescopes at continental baselines from CHIME (Leung et al. 2021). 
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APPENDIX 


A. CANDIDATE REPEATER SOURCES 


As we have established that narrowband emission and downward-drifting subbursts are ubiquitous among repeater 
bursts, we select candidate repeater sources among the sample of one-off events in Catalog 1 that also show these 
characteristics and that have had relatively low exposure. These sources can be good candidates for follow-up observa- 
tions. To deem an event a candidate repeater source, we require the subburst drift rate Av. < 30 MHz ms~! (in turn, 
this requires the number of subbursts to be > 2), the bandwidth of each subburst < 300 MHz and the exposure to 
the source in CHIME/FRB’s upper transit < 40 hr (which allows for the burst rate > 0.05 hr^! — roughly the lowest 
observed rate for a CHIME/FRB repeater; Fonseca et al. 2020). We find eight candidates based on these criteria. 
Although in Catalog 1 we only analyze the highest S/N beam detection of each event, we inspected the intensity data 
for the surrounding beams for each candidate repeater and we exclude two candidates for which we see evidence for 
more broadband emission by comparing multiple beams (see 82.1.3). In Figure 8 we compare subburst separation and 
CHIME/FRB source exposure, as well as subburst bandwidths, for pairs of subbursts in Catalog 1. Table 2 lists the 
six one-off events in Catalog 1 that are potential repeater sources based on the aforementioned criteria. 


€ 4 One-off events € o Broadband sub-burst pair © Repeater candidate 
O O Repeater bursts € o Narrowband sub-burst pair 


10! 102 
Upper transit exposure [hr] 


Figure 8. Subburst pair separation in frequency and source exposure (in the CHIME/FRB upper transit). Pairs of subbursts 
from one-off events (green/gray filled diamonds) and repeater bursts (orange/gray open circles) are shown. Multiple pairs of 
subbursts associated with the same event are linked with verticallines. Non-gray markers indicate narrowband pairs (bandwidth 
< 300MHz). A source with a Poisson burst rate of 0.05 hr ! (gray dotted line) — roughly the lowest observed rate for a 
CHIME/FRB repeater (Fonseca et al. 2020) — is expected to burst at least twice (i.e., is observed to repeat) in 40 hours. A 
one-off event with > 40 hours of exposure would thus have been detected as a repeater if the source would be as prolific as the 
known repeaters. Repeater candidates are identified with open purple diamonds. 


18 


Table 2. One-off events in Catalog 1 that are good candidates for follow-up observations to detect repeat bursts. Parameters 
with their uncertainties (68%; in parentheses) are from Catalog 1. 


FRB R.A. Decl. DM Exposure 
P] [] [pe cm 7] [hr] 
20181125A  147.9(2) 33.9(4)  272.19(2) 9(9) 
20190308C  188.36(3)  44.4(2) 500.52(3) 19(7) 
20190527A . 12.5(2) 7.99(7) 584.6(2) 9(7) 
20190422A  48.6(2)  35.2(2)  452.30(1) 15(8) 
20190423A  179.7(2)  55.3(2) 242.6456(6) 24(9) 
20190601C . 88.5(2)  28.47(6)  424.066(6) 11(7) 


